Fulde-Ferrell-Larkin-Ovchinnikov critical polarization in ID fermionic optical lattices 
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We deduce an expression for the critical polarization Pc below which the FFLO-state emerges 
in one-dimensional lattices with spin-imbalanced populations. We provide and explore the phase 
diagram of unconfined chains as a function of polarization, interaction and particle density. For 
harmonically confined systems we supply a quantitative mapping which allows to apply our phase 
diagram also for confined chains. We find analytically, and confirm numerically, that the upper 
bound for the critical polarization is universal: p™ ax = 1/3 for any density, interaction and con- 
finement strength. 
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Superconductivity, which has celebrated 100 years [![, 
has fascinated scientists since it was discovered Q. In 
non-polarized, fermionic systems, superfluidity is de- 
scribed by Bardeen-Cooper-Schrieffer (BCS) theory Q, 
in which attractive fermions with zero-momentum are 
paired. External magnetic fields, as well as polarization 
induced by spin-imbalanced populations, are expected to 
destroy the BCS-pairing mechanism and force the system 
into a polarized normal Fermi liquid. However, according 
to Fulde and Ferrell Q and, independently, to Larkin and 
Ovchinnikov Q, there would exist a polarization regime 
for which superfluidity still survives against the normal 
regime, referred to as FFLO phase. This exotic coexis- 
tence of superconductivity and magnetism is predicted to 
manifest by a spontaneous breaking of spatial symmetry, 
characterized by an oscillating order parameter [g| . 

Experimentally, there are only indirect evidences of 
the FFLO phase, in solid-state materials Q and in ID 
tubes [sj. Theoretically, considerable progress has been 
achieved in understanding the FFLO general properties 
asm , regarding the complexity of dealing with many- 
particle interactions, large systems, and the harmonic 
confinement necessary to describe state-of-the-art exper- 
iments. The regime of polarizations at which the FFLO- 
phase can be found has been determined however only on 
empirical grounds: for specific systems and parameters, 
by means of FFLO witnesses. 

We here derive an expression from general, system- 
independent considerations for the critical polarization 
Pc below which the FFLO-state emerges. We depict 
the phase diagram for unconfined systems and provide 
a quantitative mapping which allows its application also 
to harmonically confined chains. Although the critical 
polarization exhibits a strong parameter dependence, we 
find that its upper bound is universal: P™ ax = 1/3, for 
any density, interaction and confinement strength. 

Our method consists in solving the one-dimensional 
Hubbard model (25j . 

H = -t c\ a c ja + U^n^hii + V^xfniv, (1) 

<ij>,cr i ia 

at zero temperature, with onsite attractive interaction 



U , next-neighbor tunneling t, harmonic confinement of 
strength V, lattice size L, total number of particles 
N = N t + Ni, with density n = N/L = n t + n i 
and spin-imbalance P = (Nf — N^)/N. Furthermore, 
c\ a (cia) is the fermionic creation (annihilation) oper- 
ator, fi^ — c\ a Ci a counts the particle number on site 
i, and a runs over the spin values t and J,. Ground- 
state properties of finite chains are obtained via Density- 
Matrix Renormalization GroupjDMRG) [27| and Den- 
sity Functional Theory (DFT) [28| calculations, the latter 
with a local spin-density approximation to the exchange- 
correlation energy |29J. A Bethe-Ansatz (BA) solution 
is used for average properties of infinite chains, through 
a recently derived analytical parametrization (BA-FVC) 
26]. 

Let us first consider the spin state occupation proba- 
bilities. Each chain site has a four-dimensional Hilbert 
space: paired spins, single occupation of a spin-up or 
spin-down particle, and vacuum. The key quantities are 
the paired probability, 



1 dE {n,P,U) 



dU 



(2) 



closely related to superfluidity, and the majority un- 
paired probability, 
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) ~ W U = 77 [ 1 + p ] - w tl- ( 3 ) 



associated with the normal phase, where Eo(n, P, U) is 
the ground-state energy. 

We start, in Figure 1, by monitoring w-fj. and w-|- as a 
function of the spin-imbalance in the strongly interact- 
ing regime (U = —8t). The probabilities are found to be 
independent on the confinement in this regime, what is 
plausible, since the harmonic confinement plays a less im- 
portant role for strong interactions. Also, as w-j-j, and w-|- 
are average quantities, we observe that the usually less 
precise methods, BA-FVC for unconfined and DFT for 
confined systems, show very good agreement with the 
more precise DMRG data. In particular, for the inter- 
section point Pj = P(w-j- = Wf4_), Table 1 shows that 
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FIG. 1. Occupation probabilities (unpaired w-f and paired 
Wfj,), averaged over the entire lattice, as a function of polar- 
ization P. The intersection point Pi coincides for different 
computational approaches, for confined (V — 0.02i) and un- 
confined chains (V = 0), in the strongly interacting regime, 
U — —8t, and at particle density n = 0.7: analytical BA-FVC 
results for the infinite chain, DFT and DMRG calculations for 
finite chains (L = 80) with open boundary conditions. 



this good agreement holds for any density and interac- 
tion strength. The rather good performance of BA-FVC 
and DFT allows us to obtain Pj analytically for uncon- 
fined chains (via BA-FVC) and numerically, at about l%c 
of the computing time of the DMRG calculation, for con- 
fined systems (via DFT). 

At a first glance, Fig.l suggests no special feature for 
the probabilities in the entire range of polarization: they 
seem to vary linearly with P. If this was the case, the 



TABLE I. Critical Polarization Pc for several densities and 
interaction regimes, confined and unconfined chains, obtained 
with DMRG and DFT, for chains with open boundary con- 
ditions and L = 80 sites, and by BA-FVC [26j, for infinite 
chains. 
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FIG. 2. a) Pair and majority unpaired susceptibilities and 
b) minority species susceptibility, for n = 0.7 and U = — 8t, 
extracted from the analytical BA-FVC results, which clearly 
exhibit an inflection point at P ~ Pi. This corresponds to the 
critical polarization Pc which delimits the superfluid phase, 
for P < Pc (FFLO), and the polarized normal phase, for P > 
Pc (PN). The insets show the same quantities for U = — 3t: 
as the pairs are less bound, Xpair varies in a larger range (here 
Pc{n = 0.7, U = -3t) ~ 0.14). 



pair susceptibility Xpair — \dw^/dP\ (or equivalently, 
the majority species susceptibility xt = \dwf/dP\) would 
be constant in P, what means that there would exist no 
special regime of polarization in which the pairing mech- 
anism, or equivalently superfluidity, is protected. Figure 
lb instead reveals that Xpair is not a constant: On the 
contrary, it clearly shows two distinct phases with in- 
flection point precisely at the intersection point Pj. For 
P < Pf the pair susceptibility decreases with imbalance, 
i.e., the system acts against the increase of P, to pro- 
tect the pairing mechanism — what is consistent with the 
FFLO superconducting phase predicted behavior. On 
the other hand, for P > Pj, the unprotected pairs are in- 
creasingly more susceptible to P, i.e., the system acts in 
favor of enhanced polarization — characteristic of a nor- 
mal, non-superfluid phase Similar behavior is also 
observed for weaker interactions, as shows the inset of 



Polarization 

FIG. 3. Phase diagram as a function of polarization P 
and attractive interaction with strength |f/|/t, for unconfincd 
chains (for V 7^ the demarcation line between FFLO and 
PN is shifted according to Fig.5). Yellow (light gray) and 
green (dark gray) areas represent superfluid (FFLO) and po- 
larized normal (PN) fluid, respectively. P = identifies the 
BCS phase, while P = 1 the fully polarized phase (FP). The 
demarcation line between FFLO and PN, obtained via Eq.@, 
is shown for four densities. The vertical dashed line delim- 
its the analytically found upper bound p™ ax = 1/3. Stars 
represent the situations depicted in Fig. 4, panels (a) — (d). 



ti , 1 , U o.O W , 1 . 

20 40 20 40 

Lattice site Lattice site 

FIG. 4. Evolution of magnetization (m = rpf —n±), unpaired 
(w-f ) and paired (w-f-j,) probability profiles with imbalance P: 
(a) and (b) in the FFLO regime (P < P c ); (c) and (d) for 
a polarized normal phase (P > Pc)- Solid lines indicate the 
averages of Wfj, and Wf over the entire lattice. The profiles 
are obtained via DMRG for L = 40, TV = 23, U = -8t and 
V = 0. P c (n = 0.58, U = -8t) = 0.29, and each panel is 
located — by an asterisk — in the phase diagram in Fig.3. 



Fig.lb for U = -3t. 

The pairing mechanism in the FFLO regime is pro- 
tected by the unpaired minority species, wj, (for P > 0). 
When increasing the imbalance (for constant total num- 
ber of particles), the system has two flipping channels: 
either (I) from an unpaired state, |i) — > |t), or (II) from 
a paired state, |tt) — > It) > It)- Energetically, channel (I) 
is favored, with energy cost restricted to the polariza- 
tion, while in (II) there is the additional cost of breaking 
a pair. So one would expect that channel (I) is preferable 
for any P. However, since the probability of finding the 
state \l) is very small (though finite) for all P (due to 
the attractive interaction), spin flips always have contri- 
butions from both channels (I) and (II). Note that chan- 
nel (I) does not affect pairing, while channel (II) does. 
Consequently, the larger the weight of (I), determined by 
the susceptibility \i = \dwi/dP\, the more robust are 
the pairs. This is supported by Fig. 2b: For P < Pj, 
Xx increases with polarization, while for P > Pi it de- 
creases with P increasing, due to the minority species' 
then almost negligible population. Beyond Pc, polariza- 
tion is therefore progressively created by the breaking of 
pairs, channel (II), and the superfluid component must 
fade away. 

Thus Pj has an important physical meaning: it cor- 
responds to the critical polarization Pc below which the 
FFLO-state emerges, and is given by 



P c (n,U)=± 



4w n (n,P c ,U) 



(4) 



where we have set w-m = w-j- in equations ([I]) — ©- The 
sign in Eq.Q depends on the majority species, spin up 
(+) or down (— ). Figure 3 shows the phase diagram as a 
function of interaction and polarization, with the demar- 
cation line between FFLO and polarized normal phase 
(PN) (defined by Eq.@) indicated for various densities 
n. We find that for small U and n the FFLO-phase is sup- 
pressed: already small, finite values of P suffice to induce 
a transition directly from the BCS to a partially polar- 
ized normal regime (ip| |. For larger densities and weak 
interactions, the FFLO area increases with n, whereas it 
is nearly independent of n for larger values of U. That 
because the larger the density the larger the initial (at 
P = 0) pairing probability, which helps the superfluidity 
to prevail for higher polarizations at moderate U, but 
becomes irrelevant for strong interactions. 

In fact, one can analytically obtain this upper bound 
for the critical polarization, P™ aa: , by calculating the 
maximum value of w^i (see Eq.(4)). By applying usual 
particle-hole transformations between attractive and re- 
pulsive systems [30, |3l| , 

TI 

w n (n, P, U < 0) = - [1 - \P\] - w n (n', P', \U\), (5) 

where n' = nP + 1 and P' = (n — l)/n', we see that 
the maximal value w^| aa; ([/ < 0) corresponds to the 
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FIG. 5. Mapping between confined and unconfined chains 
via DFT calculations: Shift of the critical polarization, 
Pc{n,U,V) — Pc{n,U,V = 0), as a function of interaction 
\U\/t and confinement strength V/t, for different particle den- 
sities. L = N = 80 and U = -St. 



minimum of w-f4.(|?7|). This is precisely the case at 
U — > — oo: double occupancy in the repulsive system 
vanishes, Wf4.(|C7| — > oo) = 0, and double occupancy in 
the attractive system, Eq. (5), reaches its maximum, 
w™ aa; (C/ -> -oo) = (ra/2)(l - |P|). Plugging this latter 
in our formula (4) for Pq, we find the upper bound 



\P™ ax \ 



1 

3' 



(6) 



which is completely consistent with the phase diagram 
(Fig. 3) and with our numerical observations for large U 
(Table 1). 

In the limit of very strong attractive interactions, U — > 
— oo, Eq.© reads 



w n (n, P,U -> -oo) = - [1 



\P\ 



(7) 



since in this limit the double occupancy for the repulsive 
system vanishes, Wfi(n', P', \U\ — > oo) = 0. Thus, the 
majority-spin probability for P > (Eq.(3)) is given by 



w^(n, P, U 



— oo) = nP 

_ N /N t -N 
~ 1 V N~ 
= rif- m, 



(8) 



where m = — is the magnetization. This conver- 
gence of the averages w^ and m in the limit U — > — oo 
suggests that, for strong but finite interactions, the gen- 
eral properties of the unpaired probability profile can be 
observed in the local magnetization. Indeed, in Figure 4, 
we verify numerically the local agreement between w^-j- 
and rrii, for U — —8t. In the superfluid regime, P < Pc 
(Figs. 4a, 4b), we find that w-^j. dominates the entire 
chain, while the unpaired spins are inhomogeneously dis- 
tributed, with an accumulation (local probability above 
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FIG. 6. a) Critical polarization Pc as a function of the 
confinement strength V, in the regime of weak interactions 
and small densities where, at V = 0, there is no FFLO-phase 
(Pc ~ 0, also see Table I and Fig. 3). As V increases, the 
FFLO-state emerges (finite Pc)- Inset: The effective density 
at the potential core, as defined in Ref.[l5|. in units of the 
average density n, as a function of V. b) Critical polarization 
of confined chains, evaluated: i) for the entire chain (P at 
Wn = w t ), ii) at the core (P core at wff e = w| ore ) and iii) 
for the entire chain, though at the onset of the normal phase 
in the core (P at w^° re = w' 
analytical upper bound P ( 
U = -81 



corey rj^e solid line indicates the 
g ax = 1/3. Here L = N = 80 and 



the average) at the center of the chain. In the normal 
regime, P > Pc (Figs. 4c, 4d), the opposite occurs: w-j- 
prevails over the pairs in the entire chain, with Friedel- 
like oscillations around its average, while w-fj, accumu- 
lates at the chain center. 

Experimentally, ultracold atoms are among the most 
suitable systems for the detection of exotic superfluid- 
ity [1, H, El, IHl- For strongly interacting lattice sys- 
tems (\U\ > 8t), we show, however, in Fig.l and Table 
1, that the experimental (22| harmonic confinement V 
only slightly shifts the critical polarization. Therefore, 
the phase diagram in Fig. 3, constructed for unconfined 
chains, can be directly applied for harmonically confined 
systems in this regime of interaction. For weakly inter- 
acting systems, though, depending on the actual particle 
density, V implies a significant shift of Pc which cannot 
be neglected anymore. 

Hence, by comparing the critical polarizations of con- 
fined and unconfined chains, obtained with DFT (similar 
to Fig.l, but now for several n and U), we provide, in 
Figure 5, the necessary shift on Pc, which allows one to 
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apply our phase diagram for any harmonically confined 
system. Independently of the required shift, we find that 
the critical polarization is still limited to P™ ax = 1/3 for 
any density, interaction and confinement strength, i.e., 
the analytical upper bound obtained here is universal. 

Finally, we remark that this upper bound, which is con- 
sistent with the original prediction and also with quan- 
tum Monte Carlo calculations [3], has an apparent dis- 
crepancy with recent investigations in harmonically con- 
fined systems pH - [l7| . which report FFLO at higher po- 
larizations, as much as Pc ~ 0.8. This is however a sim- 
ple artifact of the confinement: While unconfined chains 
present no FFLO-state (Pc = 0) for small n and U, 
Figure 6a reveals that, depending on the intensity V, 
harmonically confined systems, conversely, do have an 
FFLO phase (P c ^ 0) (see also DMRG data for U=-2t 
in Table I). That because V reduces the volume where 
both species can be found (the central core) [HI, , such 
that the core density actually is much larger than the 



average n (see inset of Fig. 6a) and thus the FFLO-state 
emerges. This smaller effective volume is also responsi- 
ble for the observations of FFLO at polarizations beyond 
P™ ax pl| - [l7j : while the FFLO witnesses were observed 
exclusively at the central core, they were associated with 
a larger global polarization (defined for the entire chain) , 
resulting in values of Pc much larger than our upper 
bound p™ ax = 1/3. This apparent discrepancy disap- 
pears though whenever one compares the FFLO signa- 
tures (in our case w|° re and w| ore ) to the polarization 
within the core, as we do in Figure 6b. 
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